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All references are textbooks and published papers 
Reference(s) used listed on each slide, lower left, as ref# 








NASA Themes and Libration Orbits 

Enterprises of Space Sciences (SSE) and Earth Sciences 
n ar e a combination of several programs and themes 



• The Triana mission is the lone ESE mission not orbiting the Earth. 

• A major challenge is formation flying components of Constellation-X, MAXIM, TPE, and 

Stellar Imager,_ _ - 3 — 




















Earth Science Launches 
Low Earth Orbit Formations 



- Parasol (04) 

- OCO (tbd) 
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SPECS (Submillimeter Probe of the Evolution of Cosmic Structure): 
Interferometer 1 km at L2 




Future Mission Challenges 
Considering science and operations 



> Solar sail applications 

> Continuous control to reference trajectories 

> Tethered missions 

> Human exploration 









Background on perturbation 
theory / accelerations 
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FORCES ON 
PROPAGATED ORBIT 



# 3-7 



Gaussian Lagrange Planetary Equations 
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Potential Accelerations 



-3J 2 nR e 2 r. 

















Contains lag heating terms 
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Solar Radiation Pressure Acceleration 











Other Perturbations 
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#3-7 



oo 



Inertial pointing spacecraft could have drastic changes in B c over an orbit 


Ballistic Effects 

Varying the mass to area yields different decay rates 

Sample: 100kg with area of 1, 10, 25, and 50m 2 , C d =2.2 



satlOOtoI.EIapsedDays (Days) 



Solutions to ordinary differential equations (ODEs) to 
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Matlab ODE 4/5 is a variable step RK 
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ECEF - Navigation 
Topographic - Ground station 






Describing Motion Near a Known Orbit 


CD 

CD 



CD 

c/3 

O 

h-> 

CD 

O 

Oh 

P 

C/3 

C/3 

• T“H 

i 5—i 
CO 

<4-H 

• i-H 

Vh 

<D 

> 


C/3 

<D 

CD 

C/3 

C/3 


CD 


£ 

o 


+ 


l 1?^. 

co> CO 


5-< 

o 


K & 

• H 

4 —* 

O cd 
cd *r< 
X ^ 

CD 

73 (D 

•T-( 1-^J 


+ 

* 

Ik. 

II 

+ 

* 

Ik. 

II 

£ 

© 


C/3 

w 

o 

o 

C-M 

o 

4-* 

<D 

C/3 

?h 

cd 

CD 

P 

• t-H 
i H 

cd 

73 

2' 

T3 

p 


I 

% 


I'w I*. 

co^l CO 
+ 

* 

Ik 




lf<. 

I 

w 


t—H 

cd 

CN 

S -^ 

PH 

bX) 

C/3 

X 

P 

Td 

p 

cd 

-4—» 

<D 

(D 

X) 

• rH 

4-> 

P 

H—* 

-J--* 

C/3 

£ 

73 

4—> 

• t-H 

cd 

-9 

C/3 
<1 1 

© 

1 C+H 

P 

m 


^3 



This is important since it will be our starting point for everything 
that follows 




- Hill’s - RTN 

- Clohessy-Wiltshire - RAC 

- LYLH - RIC 







Describing Motion Near a Known Orbit 
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(N 



# 5,6,7 




Transforming the EOM 













Transforming the EOM yields 
Clohessy-Wiltshire Equations 




a . (associated with T) has twice 

cos(nt) + — sin (nt) the amplitude of the x motion 



Relative Motion 

A numerical simulation using RK8/9 and point mass 
Effect of Velocity (1 m/s) or Position(l m) Difference 



- 0.045 - 0.040 - 0.035 - 0.030 - 0.025 - 0.020 - 0.015 - 0,010 - 0.005 0.000 

mairvAlongArcTrackSeparatlon (Km) 







Graphics Ref: Collins, Meissinger, and Bell, Small Orbit Transfer Vehicle (OTV) for On-Orbit 
Satellite Servicing and Resupply, 15th USU Small Satellite Conference, 2001 












Consider two trajectories r(t) and R(t). 



# 5 , 6,7 




Lambert Problem 
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‘single pass’. 





















EO-1 GSFC Formation Flying 

New Millennium Requirements 



re re 



The onboard flight software, called the EFF, shall provide the interface between the ACS / 
C&DH and the AutoCon™ system for Autonomy for transfer of all data and tables. 














Observation Ov! 









EO-1 Formation Flying Algorithm 
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where F and G are found from Gauss problem and the f _& 9 
series and C found through universal variable formulation 
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Inclination Maneuver Validation : Computed AV at node crossing, of 
~ 24 cm/s (114 sec duration), Ground validation gave same results 






Difference in EO-1 Onboard 
and Ground Maneuver Three-Axis AVs 



# 10,11 
























Formation Data from Definitive Navigation 

Solutions 



formation maneuvers 
maintained to 3km 








Formation Data from Definitive Navigation 

Solutions 
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Formation Data from Definitive Navigation 

Solutions 











EO-1 Summary / Conclusions 
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Summary / Conclusions 
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Formation Flying and Multiple Spacecraft Missions 
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NASA Libration Missions 
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• MAXIM Lissajous, Formation Flying of Multiple S/C 

• TPF Lissajous, Formation Flying of Multiple S/C 

(Previous missions marked in blue) 












# 1,12,13 Farquhar et al [1985] Trajectories and Orbital Maneuvers for the ISHE-3/1CE, Comet Mission, JAS 33, No. 
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Spacecraft: Mass=818kg, Three Axis Stabilized, 

Notable: First Gravity Assisted Constrained L2 Lissajous Orbit; Map-earth Vector 

Remains Between 0.5° and 10° off the Sun-earth Vector to Satisfy 
Communications Requirements While Avoiding Eclipses 52 
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Ecliptic Plane 
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Aspherical Configurations (Position & Rates) 





















CR3BP Analysis of Phase Space Eigenstructure 

Near Halo Orbit 



Solution to Variational Eqn. in terms of Floquet Modes: 

8* (0 = (/) = Y J c j {t)e j (f ) = E(t)c 






Natural Formations: 

Quasi-Periodic Relative Orbits —> 2-D Torus 









Floquet Controller 

(Remove Unstable + 2 of the 4 Center Modes) 
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Deployment into Torus 
(Remove Modes 1, 5, and 6) 
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Deployment into Natural Orbits 
(Remove Modes 1, 3, and 4) 
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Natural Formations: 
Nearly Periodic Relative Motion 
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Natural Formations: 
Expanding Vertical Orbits 
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Fixed Relative Distance and Orientation 
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Formations Fixed in the Inertial Frame 
























Formations Fixed in the 
Rotating Frame 



y (10° km) 
















2-S/C Formation Model 
le Sun-Earth-Moon System 



Ephemeris System = Sun+Earth+Moon Ephemeris 










Goddard Space Flight Center 
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Linear Targeter 









































Achievable Accuracy 
via Targeter Scheme 



Formation Distance (meters) 






Continuous Control: 

LQR vs. Input Feedback Linearization 







r- 



n 14,15 




—» Substitute x° (/) into EOMs and solve for u d ° (V) 
-*u d (t) = u d °{t)-K(t)(x(t)-x°{t)) 

-> Apply the computed control input to the perturbed EOMs 



CO 

CO 

CD 

U 

O 

L_ 

Q. 


LL. 


M 



CO 


o 


• y—4 


1/3 


*a 


CD 


J * 

O 


cd 


5-i 


aS 

•ls<. 

43 

*\ 

o 

|Jk 

CD 

s — ^ 

C/3 

c 

&0 

o 

II 

P 

^ - v 

C/3 

o 

CD 

5-i 

K 

T3 

1 

CD 

IJk 

5-i 

* 


(Z) 

(N 

CD 

CO 

n 

CD 

1 

43 



O 

c*\ 

*|fc* 


1 

o 

• r-H 



3 s 

13 

<N 

i 

c 

aS 

1 

CD 

:hk 

a 

II 


CD 

:K 

Q 

t 

H 

1 

P 


CD 

i * 


GO 




aS 


C/3 


C/3 


CD 


CD 


CD 


33 


4-J 


33 


P 


33 


• *-H 


r—H 


O 


£ 


33 


O 

43 

(D 

CD 

43 

c3 

P 

4 —> 

'U 

4^ 

CD 

•e 

ft 

& 

•P 

"O 


33 

P 

as 

p 

CD 

o 

+-» 

33 

33 

P 

O 

• r-H 

& 

aS 

5—( 

O 

CD 

30 

CN 

CD 

CD 

+-> 

33 

E 

* rH 

4—> 

13 

33 

CD 

• i-H 

*c 

I 

+-> 

33 

O 

33 

P 

33 

43 

33 

o 

aS 

30 

CD 

CD 

4-> 

PQ 

<; 

CN 

rn 

P 

P 

0 ) 

—«■■ * 

CD 

GO 

LTj 


\A 

o 


C/> 

5-1 

<D 

'g 

03 

J-i 

OS 

43 

o 

CD 

co 

33 

O 

a, 

C/3 

CD 

5-1 


T3 

CD 

5-i 

• 

C/3 


CD 


X3 


43 


CD 

> 

<D 


43 


o 

aS 

O 



00 

r-~ 



SI # 



t-1 


C/2 

cu 


a 


C/2 


c 

H 

y—i 

| 

o 



C/5 

C/2 

'l_ 


& 

£ 

05 

o 

Q_ 

O 


E 

G\ 


o 

II 

g 

u 


i2 

—1 

CN 

IT2 

LJL 


cn 

i — i 

q 


CO 

> 

o 

o 

£ 

ttL 

o 

i/n 

i 

cy 

to 

1 1 

i 




^ E 

IK « 

Q 
























Output Feedback Linearization 
(Radial Distance Control) 



Scalar Nonlinear Constraint on Control Inputs 

h(7(t),¥(t)^-u(t) T r(t) = 0 




Output Feedback Linearization (OFL) 
(Radial Distance Control in the Eohemeris Model ') 


00 




Critically damped output response achieved in all cases 
Total AV can vary significantly for these four controllers 






OFL Control of Spherical Formations 
in the Ephemeris Model 



Relative Dynamics as Observed in the Inertial Frame 



















OFL Controlled Response of Deputy S/C 
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OFL Control of Spherical Formations 
Radial Dist. + Rotation Rate 
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Inertially Fixed Formations 
in the Ephemeris Model 



inertially fixed formation pointing vector (focal line) 






Nominal Formation Keeping Cost 
(Configurations Fixed in the RLP Frame) 
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Max./Min. Cost Formations 
(Configurations Fixed in the RLP Frame) 




Nominal Relative Dynamics in the Synodic Rotating Frame 




Formation Keeping Cost Variation 
Along the SEM and L 2 Halo Families 
(Configurations Fixed in the RLP Frame) 



x (10 km) A (10 6 km) 



















Continuous Control in the Ephemeris Model: 
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Nearly periodic & quasi-periodic formations in the RLP frame 
Floquet controller: numerically ID solutions + stable manifolds 





Some Examples from Simulations 
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Periodic Reference Orbit 
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Centralized LQR Design 



#16,17,20 








Centralized LQR Design 











Disturbance Accommodation Model 
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Disturbance Accommodation Model 
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Motion of Formation Flyer With Respect to 
Reference Spacecraft, in Local (S/C-l) Coordinates 
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Stellar Imager Concept 

(Using conceptual distances and control requirements 
to analyze formation possibilities) 
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spacecraft, and environmental conditions 

#16,17,20 
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Stellar Imager Mission Study Example Requirements 
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Origin in figure is Earth 
Solar rotating coordinates 























Algebraic Riccatic Eq. *4 S(BR~‘B t )S -SA-A t S 
time invariant system 
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Simulation Matrix Initial Values 








Results - Libration Orbit Maintenance 
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Results - Libration Orbit Maintenance 
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Represents both 0.5 and 4 km focal lengths 




























Results - Formation Slewing 
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Hub estimation 0.5 km Drone estimation 0.5 km 

separation / 90 deg slew separation / 90 deg slew 































6.3662 0.9577 0.9773 

7.1135 1.1406 1.1534 

10.4112 2.2331 2.2357 





Formation Slewing Average AVs (without noise) 



0.8864 0.2192 0.2360 

2.4781 0.6010 0.6566 

7.8182 1.9822 1.9846 




Results - Formation Reorientation 
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Results - Formation Reorientation 
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Results - Formation Reorientation 
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Summary 

(using example requirements and constraints) 
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control; ior m couia be accomplished with better sensors to 
lessen the effect of the process and measurement noise. 



Summary 

(using example requirements and constraints) 
the controller and varying the maneuver intervals 
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• Nonlinear control and estimation 

• Second and third control tiers and new control 
strategies and algorithms 
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DRO Mission Metrics 
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Distant Retrograde Orbit (DRO) 



Asymptotic Orbit in 
the Restricted Three- 
Body Problem 



Earth Distant Retrograde Orbit (DRO) Orbit 



S1RA_DR0 






r# 1,18 



















































DRO Formation Control Analysis 



# 1,18 





DRO Formation Control Analysis 
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DRO Formation Control Analysis 



#1,18 




Formation Control Analysis 



d) Resize from 25km to 50km 

e) Maintain 50km sphere (strict PD control) one month 

f) Maintain 50km sphere (loose control) one month 
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DRO Formation Control Analysis 

> Earth/Moon L4 Libration Orbit 

> Spacecraft controlled to maintain only relative separations 
^ Plots show formation position and drift (sphere represent 2 



0.000 0.025 0.050 0.075 0.100 0.125 0.150 0.175 

r# 1,18 _ AI.EIapsedDays (Days) _ 






General Theory of Decentralized Control 
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DST/Numerical Comparisons 





ymplectic Integration Schemes? 



Libration Point Trajectory Generation Process 























































MAP Mission Design: DST Perspective 



irajectory ueneratec 
with Manifold States 
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Coordinate system transformation 
Orbital to inertial coordinates 








Principles Behind Decentralized Control 






The LQG Decentralized Controller Overview 
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TRANSMISSION VECTOR 

TO NETWORK 





